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Abstract 

Presented here is a matrix inversion method utilizing quantum searching algorithm. In this 
method, huge Hilbert space as a whole spanned by myriad of eigen states is searched and 
evaluated efficiently by sequential reduction in dimension one by one. Total iteration steps 
required for search are proportional to the number of unknown variables. Our method could 
solve very large linear equations with sufficiently high probability faster than any existing 
classical algorithms, which roughly depends on the cube of unknown variables. 
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1 Introduction 



2 Quantum matrix inversion 



Since GroverQ] has first invented his quantum 
searching algorithm quadratically faster than clas- 
sical counterparts, several applications of the al- 
gorithm have been considered, for example, find- 
ing the minimum[2|, estimating the median[P, 
quantum counting|4], collision problem]^, undi- 
rected graph conncctivity|6| and protein sequence 
comparison^. On the other hand, this searching 
algorithm is thought to be unsuitable for such mas- 
sively number crunching computer simulations like 
aerodynamics, car crush, electric circuitry, weather 
prediction, molecular biology, finance and so on. 
Frequently, large part of those simulations are de- 
voted to solve linear equations, which are mathe- 
matically written in matrix form. In this paper, 
we present a method to solve matrix equation effi- 
ciently on quantum computer by utilizing quantum 
searching algorithm. In our method, huge Hilbert 
space as a whole spanned by myriad of eigen states 
is searched and evaluated efficiently by sequential 
reduction in dimension one by one. 

In the next section, we describe the matrix in- 
version method both naive and improved in detail. 
The third section estimates the expected computa- 
tional performance numerically. The final section 
concludes this paper. 



Any system of linear equations can be written as 
single matrix equation as following 



Ax 



0. 



where A denotes matrix with n-columns and m- 
rows, and h denotes vector with n-rows. For 
simplicity, we assume that the system is neither 
overdeterminate nor underdeterminate. This as- 
sumption imposes n = m and det(^) on ma- 
trix A and guarantees unique solution. The prob- 
lem to be solved is to locate x satisfying above 
equation for given A and b. This kind of prob- 
lem has been investigated for many years and the 
baseline algorithm was discovered by great Ger- 
man mathematician, Carl Friedrich Gauss, in nine- 
teenth century. Unfortunately, his algorithm some- 
times falls in instability and inaccuracy for a large 
matrix on digital computers mainly because of ac- 
cumulation and amplification of truncation error. 
Although many improvements, enhancements and 
innovation have been persued on classical algo- 
rithm, comutational workload propotional to the 
cube of unknown variables is still required. On the 
other hand, quantum algorithm utilizes novel char- 
acteristics of quantum mechanics known as quan- 
tum parallelism to accelerate its computation dra- 
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matically. Especially, Grover's||l[ quantum search- 
ing algorithm is quite versatile to use and already 
applied in several problems. In the following sub- 
sections, we propose the matrix inversion method 
based on his algorithm to obtain further accelera- 
tion from classical counterparts. 

2.1 Naive implementation 

In this subsection, we describe the straightforward 
implementation of Grover's algorithm on matrix 
inversion. For the first time, we prepare 3n null 
registers |0) 

10) = 



Operating Walsh-Hadamard gates on the first 
n registers provides uniformly superposed initial 
state \ip{0)), which forms discrete n-dimensional 
searching space as following 



1^(0)) 



M 



0) 



n,fn ^ 



M 

E 



where H and / denote Walsh-Hadamard gate and 
identity operator, and M denotes the number of 
numerical points in each dimension. Therefore 
IV'(O)) forms M" dimensional Hilbert space. Suc- 
ceedingly, we multiply the first row of matrix A 
with X. According to Vedral et. al. there exist 
unitary operators Un, U12, • • • , Uin for each matrix 
elements an, ai2 ■ ■ ■ ai„ which map xi, X2 ■ ■ ■ Xn to 
aiiXi, ai2X2 • • • ainXn such that 

Urj\xj)\0) = \x.i)\a,jXj) . 

After exertion of these operators, we obtain 

• ®;7i„®/«"|V(o)) 



IV'(i)) = c/ii®c/i2 
1 




Ifllria^ri) |0) 



where we introduced an abbreviation x, which de- 
notes a set consists of all n dimensional numerical 
points, that is, x = {xi, X2, ■ ■ ■ , Xn : 1 < a;^ < 
M for each i}. Then, we take summation of the 
middle registers to obtain /i defined as 

n 

/i = E "'^3^3 + ^1 ' 

where bi denotes the first element of vector b. 
Since addition also can be realized by some uni- 
tary operators B, we can calculate above equation 



through quantum gates. Here we obtain following 
state 




■•|aina;„)|/i)|0)---|0) . 

Obviously, the middle registers are useless for the 
rest of calculation, so that we clear up such garbage 
and recycle them in the next calculation step. 
Garbage erasure is realized by backward opera- 
tion, which has first devised in connection with 
reversible computer by Bennett After garbage 
erasure we have 



\m) = 




|0)|/i)|0)---|0) 



Repeating above operations for the rest of n 
rows, eventually we obtain 



IV'(3n)) = 



This finishes preparation. From here, Grover's 
quantum searching algorithm is invoked. To utilize 
his algorithm, we must identify oracle C which de- 
termine whether the argument satisfies given con- 
straints. Oracle must have a property like that 




C 



otherwise 







However, particular construction of above oracle is 
left for future work currently. Meanwhile, Grover 
showed that only unique solution survives after a 
number of iterations of state rotation. Necessary 
iteration count has estimated exactly by Boyer et 
al. (Hereafter we refer it as BBHT)|1(][| as close 
to jVN, where N denotes the total number of 
candidates in search. They also found the fail- 
ure probability as In case of naive imple- 
mentation described here, N equals Af". Conse- 
quently, computational steps required for search 
amount to roughly M^. While failure probability 
is sufficiently low, the exponential dependency on 
n forces computational steps to explode, and there- 
fore makes naive implementation definitely imprac- 
tical for especially large-scale matrices. 

2.2 Dimensional reduction 

As you see in the previous subsection, searching 
such a huge Hilbert space as a whole would be a 
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desperate effort even if quantum computer were 
available. In this subsection, we employ "a box in a 
box" strategy which, in practice, divide the search- 
ing space into a sequence of lower dimensional sub- 
spaces one by one. Basically procedures required 
here are similar to the previous ones. However, 
we put BBHT algorithm instead of Grover's algo- 
rithm immediately after the calculation of |7/;(3)), 
since BBHT can find out multiple solutions. In 
each iteration step, we use oracle C such as 

C'^l ^ if /i = 

1^ otherwise ' 

instead of C. According to BBHT, multiple so- 
lutions are obtained after approximately j^^N/t 
iterations with failure probability t/N, where t 
denotes the number of solutions. Provided that 
the following inequality similar to diagonally dom- 
inance holds for any i 

2max|a,j| < ^ |ay| + , 

3 

then t always forms n—1 dimensional complete, 
i.e. undipped, intersection. In such cases, we can 
assure that TV = M" and i = M"-i. Therefore, 
iteration count and failure probability in this stage 
can be estimated as approximately jVm and 1 /M 
respectively. After the search in the first dimen- 
sion, we obtain the following state 

x' 

|0)...|0)|/i)|0)...|0), 

^ V ' ^ V ' 

n n — 1 

where x' denotes a subset of x, that is, x' — 
{xi,X2, • • • , : /i = ; 1 < Xi < Af for each i}. 
In other words, only n—1 dimensional subspace 
survives through the first searching stage. Next, we 
multiply the second row of matrix A with surviv- 
ing superposed state x' as just like the way we did 
previously. Succeedingly, summation and garbage 
erasure follow. Then we put BBHT again. Oracle 
to be used here is 

^„^f 1 if/2 = 

1^ otherwise ' 

instead of C". This time, volume of searching space 
is M"^^ and number of solutions are M"~^. Thus, 
necessary iteration count and failure probability 
are the same as before, that is, jVM and 1/M 
respectively. After the search in the second dimen- 
sion, the state is 

W5)) = ^EM---M 

x" 



|o^^l/l)l/2)|o^Ho), 

n n—2 

where x" denotes a subset of a;', that is, x' — 
{xi,X2, - ■ ■ ,Xn : /i = 0; /2 = 0; 1 < Xi < 
M for each i\. After repetition of this procedure 
n — 2 times for the rest rows of matrix A, search- 
ing subspace reduces to zero dimension. In this 
way, finally we get to the unique solution 

Hf^nal) = '""°«) • • • 

|O^HO)|/l)---|/n). 

n 

Here, we measure the registers to read out solu- 
tions. We can readily see that the summation of 
all iteration count amounts to ^n^/M. To get 
to the true solution, all of n searches should be 
succeeded. Thus, success probability for the total 
procedure equals to multiplication of success prob- 
ability in each search, that is, (1 — 1/Af)". Ob- 
viously, iteration count scaling with n shows ex- 
ponential acceleration from naive implementation. 
This fact shows that dimensional reduction tech- 
nique described here works well for the large ma- 
trix inversion. 

3 Performance 

In this section, we will investigate performance 
issues in detail. As shown in previous section, 
number of steps required in searching algorithm 
is jUy/M. In addition to this, steps are also 
needed to calculate inner product. Regretfully, we 
don't have any idea to make them decrease for the 
time being. Meanwhile, current digital computer 
usually assigns four bytes data(j.e. 32bits), for 
each real number. Here we also adopt this as typi- 
cal qubit size of register, that is, M = 2"^^. Now we 
can readily estimate operational count required for 
entire computation as 2n(51, 471-|-n), in which the 
factor of two reflects backward operation to erase 
garbage. On the other hand, it is well known that 
classical Gaussian elimination method requires the 
order of steps From these expressions, we 
can estimate crossover size from classical to quan- 
tum algorithm as roughly n — 321. This result 
means that even relatively small matrix can be ac- 
celerated by quantum algorithm. However, note 
that the estimation here can not be taken serious 
for the time being, since operational speed per gate 
of current digital computer is considerably faster 
than that of today's infant quantum computer. 

Another important issue to be considered be- 
sides iteration count is success probability. If the 
probability to obtain true solution is poorly low, 
we should be hopelessly exhausted recalculating so 
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many times. In previous section, the success prob- 
ability of dimensional reduction technique is writ- 
ten as (1 — 1/M)". To find out the lower bound of 
this probability, we prove following two lemmas. 

Lemma 1 Let k and n be any natural numbers 
such that k < n, then following inequality holds 

nCk < n'' . 

Proof. We prove this lemma by induction. For k = 
1, above inequality obviously holds for any natural 
number n. Provided that the inequality holds for a 
integer k such that I < k < n, LHS term for k + 1 
could be estimated as 

_ (n~k) (»- fc) ^fc . k+i 

n^k+l — , , n^k S — —n < n 



1 



1 



Notice that < n holds for any positive n and 

k. Consequently, induction principle guarantees 
above inequality for any k. ^ 

Lemma 2 Let p he a real number such that < 
p < 1 and n be a any natural number such that 
np < 1, then following inequality holds 



(1-P)"> 



1 — 2np 



I — np 

Proof. Using binary term expansion, 



n 



fc=i 



Substituting the result of lemma ^ the estimation 
continues like 



> 



k=l 



= I - 



np — (np) 
1 — np 



n-l-l 



> 



1 - 2np 



1 ~ np 
Here completes the proof. 

Substituting p = 1 /M in lemma ^ gives 



^ MJ - 1 - n/M 



Now, we estimate this lower bound of success prob- 
ability numerically. Substituting M = 2^^ as typ- 
ical number described before, we can make sure 
that success probability is no less than 93% for 



n < 2^* ~ 10^ and that the smaller n, the higher 
success probability for fixed M. 

In this section, it is shown that quantum ma- 
trix inversion enhanced by dimensional reduction 
technique could solve matrix equation with suffi- 
ciently high probability within linear time, even if 
its problem size were either relatively small or con- 
siderably large. 

4 Conclusion 

In this paper, we proposed a matrix inversion 
method utilizing quantum searching algorithm. 
The heart of our method is dimensional reduc- 
tion techneque, which is introduced to accelerate 
searching procedure. This technique enables us 
to obtain exponential speed-up over naive imple- 
mentation and sufficiently high probability of suc- 
cess. Required number of iteration steps linearly 
depends on the problem size, while quadratic steps 
are necessary for arithmetic calculations. Proposed 
method achieves acceleration over classical algo- 
rithms by an order of n. These results might open 
a possibility of quantum computer in future indus- 
trial use. Actually, our proposal is nothing but a 
primary desktop calculation and we need further 
investigation. For example, to identify particular 
gate configuration, to estimate performance strin- 
gently through emulation of gate operation and to 
understand susceptibility to error caused by trun- 
cation and decoherence are left as future works. 
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